1/i 


AD-A158  872  INVESTIGATIONS  OF  LASER  INITIATED  DISCHARGE  CHANNELS 
<U>  MICHIGAN  UNI V  ANN  ARBOR  R  M  GILGENBACH  ET  AL. 

19  AUG  85  TR-114  N00014-81-K-0700 


UNCLASSIFIED 


F/G  20/9 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-I963-A 


UNCLASSIFIED 


SECURITY  Ck«*S'*'CATIO*  O*  TMiS  PAOfi 


REPORT  DOCUMENTATION  PAGE 


tk.  RESTRICTIVE  MARKINGS 


«.  PERPORMING  ORGANIZATION  REPORT  NUMBER'S) 

Technical  Report  #114  ' 


1  MONITORING  ORGANIZATION  REPORT  NUMBS  R(E) 


.  NAME  OP  PERPORMlNG  ORGANIZATION 

University  of  Michigan 


A  OP  PICE  SYMBOL 
IlftppUnMt) 


Ta  NAME  OP  MONITORING  ORGANIZATION 

Physics  Division 
Office  of  Naval  Research 


Be.  AOORESS  (City.  S MM  mtd  ZIP  C atf»J 

University  of  Michigan 
Ann  Arbor,  MI  48109 


7A  AOORESS  (CUy.  SON  mtf  ZIP  C«B#J 

800  N.  Quincy  St. 
Arlington,  VA  22217 


Ba.  name  op  PUNOINQ/SPONSORlNG 
ORGANIZATION 


Office  of  Naval  Research 


S.  PROCUREMENT  INSTRUMENT  IOENTIPICATION  NUMBER 


Sc.  AOORESS  (City.  SUN  m4  ZIP  C«*| 

800  N.  Quincy  Street 
Arlington,  VA  22217 


IS  SOURCE  OP  PUNDINO  NOS. 


PROGRAM 
ELEMENT  NO. 


11.  Title  'ta'IvAc  Stftinly  CNmiflcatioai 

Llser  Initiated  Dis- 


1Z.  PERSONAL  autmorisi 
R.  M.  Gilgenbach 


ISA  type  op  report 

Annual  Summa 


IB.  SUPPLEMENTARY  notation 


PROJECT 

TASK 

WORK  UNIT 

NO. 

N00014 

NO. 

NO. 

81K 

0700 

i 

i 


i 


l>4 


&  T,  Kanunash 


13A  Time  covered 
pro 


r/5i 


«.  DATE  OP  REPORT  ( Yr ..  *•..  Dryl 

8/19/85 


IS.  PAO 


asawi 


COSATl  COOES 


IS  SUBJECT  TERMS  ICmNw  m  imn  if  amevy  ml  MtnMfy  By 

Channels,  laserguided  discharges 


HI.  ABSTRACT  iCmimw  an  awfw  II  awm  ad  idtmdfy  By  Mack  mumhtri 

report  progress  on  the  generation  of  intense,  microsecond, 
Qjelectron  beams.  Diode  closure  has  been  studied  for  electron  beam 
gpulselengths  up  to  4  microseconds. 

^Reduced  density  channels  produced  by  laser  induced  gas  breakdown 
jjZr^d  by  laser  guided  discharges  are  modelled  by  means  of  a  one 
«j-aimensional  gas  hydrodynamics  code  for  times  earlier  than  the 
tn set  of  convective  mixing.  Comparison  of  laser  induced  gas 

experiments  with  the  code  indicate  that  the  conversion 
C3f  laser  energy  to  gas  heating  is  limited  to  25-30%.  Col  linear^ 

(cont.) 


90  OlfT A<9wr#0*/AV Alt Afltirv  O*  AISTMCT 


91  AMTAACT  SieuAlTV  CWAM"  ICATIO* 


G  otic  users  □  UNCLASSIFIED 


SZc  OPPICS  SYMBOL 


UNClASSiPiSO'WNLIMITSO  9  same  AS  RPT  C  OTIC  USERS  Q 


Mp  name  o*  REtPONtiBLi  INOIVIDwAL 

C.  W.  Roberson 


DO  pOMM  t473,  83  ARM  edition  op  i  jan  ?»  *e  obeole te  UNCLASSIFIED 

/''Vv V->> >’*>’<* V V'j-'j. 'v‘' .•Vv'vV  .•  . 


ih  telephone  number 

(2mwr-^2 


jpsjjllE 


UNCLASSIFIED 


MCUft'TV  CLASSIFICATION  OF  THIS  FAOS 


holographic  measurements  of  channel  density  profiles  produced  by 
laser  guided  discharges  are  /found  to  be  in  good  agreement  with 
code  predictions  for  times  iess  than  50  s.  This  agrees  with 
other  experiments  which  have  indicated  onset  times  for  convective 
mixing  in  the  range  50-70  jus)  after  the  start  of  the  discharge. 


A  Monte  Carlo  Fokker  Planck  code  is  used  to  investigate  the 
scattering  loss  of  particles  from  a  relativistic  electron  beam 
propagating  in  a  low  pressure  air  channel.^  The  fraction  of 
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expression  for  force-free  propagation.  \ 
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Appendix  1 

During  the  past  year  we  have  made  significant  progress 
in  the  generation  of  multimicrosecond  electron  beams  on 
MELBA.  In  a  recent  invited  talk  at  the  IEEE  Pulse  Power 
Conference  we  reported  generation  of  4  microsecond  electron 
beams,  believed  to  be  the  longest  intense  beam  pulselengths 
in  the  0.  S.  In  these  experiments  we  discovered  significant 
effects  of  anode  conditions  upon  electron  beam  diode 
impedance.  These  results  are  fully  described  in  the  attached 
reprint  from  the  conference  proceedings. 

In  the  area  of  shock  wave  generation,  diagnosis  and 
modeling,  we  have  successfully  employed  a  1-dimensional 
hydrodynamics  code  to  estimate  the  energy  deposition  in 
channels  by  lasers  and  guided  discharges.  The  col linear 
holographic  intreferometry  results  have  been  previously 
published  and  the  comparison  of  data  to  the  code  is  given  in 
the  doctoral  dissertation  of  L.  D.  Horton,  summarized  in  the 
attached  technical  report. 

Initial  electron  beam  interaction  experiments  have 
concentrated  on  the  interaction  of  0.4  microsecond  Pebetron 
beams  with  argon  and  nitrogen.  These  results  were  presented 
at  the  1984  APS-DPP  meeting.  Some  of  the  significant 
findings  are  shown  in  Figures  1-1,  1-2,  1-3,  and  1-4.  These 
figures  give  as  a  function  of  pressure,  the  transmitted 
current  to  a  bare  collector  located  0.85  m  from  the  entrance 
foil.  Also  shown  are  the  collector  current  decay  times  for 
these  cases.  It  can  be  seen  that,  for  the  1  kA  injected 


current,  less  than  50%  is  measured  by  the  collector.  Another 
interesting  feature  was  that  the  collector  current  decay  time 
was  larger  in  argon  than  nitrogen  indicating  that  the  L/R 
time  for  argon  is  longer  (giving  a  lower  R)  for  the  postbeam 
channel . 

In  the  area  of  diagnostics  we  have  investigated  several 
new  configurations  including: 

1)  a  single  optical  fiber  probe  to  measure  spatial  and 
temporal  electron  beam  profiles  within  a  plasma,  and 

2)  Lucite  plate  diagnostics  for  measurement  of  electron  beam 
uniformity  and  potentially  emittance  in  the  diode. 

The  abstract,  from  a  paper  describing  these  diagnostics  is 
attached.  Data  from  the  fiberoptic  probe  is  presented  in  Fig. 
1-3. 

Most  of  the  effort  for  the  past  year  was  devojted  to 
improving  the  Monte  Carlo  Fokker-Planck  Code  to  incorporate  the 
self  field  of  the  beam.  This  is  in  its  last  stage  and  in  the 
near  future  the  code  will  be  tested  on  known  cases  to  assess  the 
role  of  the  self  field  on  the  scattering  losses  from  an  intense 
relativistic  electron  beam  propagating  in  a  low  pressure  laser- 
initiated  channel. 

In  the  course  of  that  study  the  code  without  the  self  field 
was  also  tested  by  investigating  the  erosion  of  an  electron  beam 
propagating  in  the  ion  focused  regime  as  generated  by  the  ETA 
experiment.  Since  there  are  several  mechanisms  responsible  for 
beam  erosion  the  usefulness  of  this  code  was  employed  to  estimate 
the  scattering  effects.  The  effects  of  multiple  scattering  were 
also  addressed  and  their  variation  with  such  parameters  as  beam 


density,  channel  neutral  and  electron  densities,  and  beam  and 
channel  radii.  The  results  of  this  study  are  summarized  in  a 
paper  published  in  the  1984  proceedings  of  the  Annual  propagation 
meeting  (copy  attached)  and  in  a  paper  presented  at  the  annual 
APS  plasma  physics  meeting. 
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MICROSECOND  ELECTRON  BEAM  DIODE  CLOSURE  EXPERIMENTS  * 

R.  M.  c> lgenbach ,  L.  D.  Horton.  R.  P.  Lucey,  Jr, 

S.  Bjdwell,  M.  Cunco,  J.  Miller,  and  L.  Smutek 

Intense  Energy  Bee*  Interaction  Laboratory 
Nuclear  Engineering  Department 
The  University  of  Michigan 
Ann  Arbor,  Michigan  4B109 


Abstract 


Experiments  have  been  conducted  on  the 
Michigan  Electron  Long  Beam  Accelerator 
1 MELBA)  to  explore  diode  closure  and  voltage 
compensation  during  collapsing  diode 
impedance.  MELBA  operates  with  parameters: 
Voltage  ■  -  0.6  to  -1.0  MVj  Current  <  60  kAj 
and  Pulselengths  exceeding  1  >is.  A  reverse 
charged  ringing  circuit  compensates  the  Marx 
generator  output  voltage  to  within  ♦  7  % 

(♦  10%)  for  pulselengths  of  1  y%  ( 1 7e  jis)  in 
a“l27  ohm  resistive  load.  Electron  beam 
generation  experiments  have  utilized  two 
types  of  field  emission  cathodes:  carbon 
brush  and  velvet  cloth.  Total  electron  beam 
pulselengths  of  4  jis  have  been  achieved. 
Voltage  compensation  is  limited  by  current 
spikes  to  ±131  (±5%)  over  the  first 
1  jis  (0.7  >ii).  Diode  impedance  and  perveance 
have  been  characterized  as  a  function  of  time 
at  several  A-K  gap  spacings.  Current 
densities.have  been  in  the  range  from  10  to 
200  A/cm  .  The  current  shows  an  initial 
step  with  a  rapid  increase  at  0.6  to  0.B  >is 
followed  by  either  current  runaway  or  a 
return  to  Child-Langmuir.  Apertured  electron 
beam  current  shows  a  similar  behavior  with 
large  amplitude  beam  fluctuations  during 
enhanced  current.  Apparent  plasma  closure 
velocities  inferred  from  diode  impedance 
data  indicate  closure  velocities  from  3.3  to 
6.2  cm/ps  dependent  upon  the  cathode  type. 

The  effective  closure  velocity,  defined  by 
the  time  required  to  short  circuit  the  A-K 
gap,  ranges  from  2.1  to  4.1  cm/ps.  Electron 
beam  uniformity  and  beam  dynamics  have  been 
investigated  by  means  of  apertured  anodes 
with  lucite  plate  diagnostics.  Voltage 
compensation  has  been  tested  by  comparing 
experimental  voltage  traces  to  transient 
circuit  computer  codes  with  empirical 
iagtedance  profiles. 

Introduction 


Most  research  concerning  intense,  pulsed 
electron  beams  has  utilized  pulselengths  in 
the  range  from  tens  to  hundreds  of  ns.  This 
was  due  in  part  to  specific  applications, 
(such  as  inertial  confinement  fusion),  which 
require  short,  intense,  pulses  of  energy. 

Accelerator  size  and  cost  can  be  reduced 
when  large  delivered  energies  are  generated 
by  Increasing  the  beam  pulselength  rather 
than  the  voltage  and  current.  Earlier 
experiments  f 1-81  have  generated  intense, 
microsecond  electron  beams  by  means  of  Marx 
generators  with  pulse  forming  lines  (PPL)  or 
pulse  forming  networks  (PPM).  In  this  article 
we  demonstrate  microsecond  electron  beam 
generation  and  voltage  compensation  with  a 
Marx  generator  coupled  to  an  Abramyan  type 
RLC  ringing,  reverse-charged  stage  (9,  10). 


The  advantages  of  the  Abramyan  circuit 
over  Marx  generators  with  PPLs  or  PPNs  are: 

1)  It  requires  fewer  circuit  elements,  thus 
providing  more  compact,  less  expensive 
generators, 

2)  Compensation  for  different  load 
impedances  can  be  changed  by  adjusting  the 
relative  charginq  voltages  rather  than  by 
changing  pfn  circuit  eleswnts, 

3)  The  Marx  generator  directly  drives  the 
load,  so  the  voltsge  requirement  is 
lowered. 

Intense  electron  beams  with  pulselengths 
exceeding  1  pi  have  a  number  of  important 
applications: 

1)  Pumping  of  lasers  in  the  visible,  uv,  and 
(possibly)  x-ray  regime  (4,10), 

2)  High  power  microwave  and  millimeter  wave 
generators  and  free  electron  lasers  til], 

3)  fusion  plasma  startup,  heating, 
and  stabilization, 

4)  Pulsed  X-ray  radiography,  and 

5)  Rapid  heating  of  materials  and  gases. 

Each  of  these  applications  imposes  different 
constraints  on  the  cathode  current  density 
and  voltage  flatness.  Electron  beam  pumped 
lasers  might  require.a  current  density  on 
the  order  of  30  A/cm  ,  whereas  high  power 
microwave  generators  operating  in  the 
collective  regime  may  need  kA/cm. 

Explosive  emission  cathodes  (1-16]  have 
generally  been  employed  in  these  applications 
because  of  their  high  current  densities  and 
immunity  to  the  hostile  plasM  chemistry  en¬ 
vironment  which  poisons  thermionic  cathodes. 
This  presents  a  problem  in  the  generation  of 
microsecond  electron  beams  because  of  the 
rapid  closure  of  the  anode-cathode  (A-K)  gap 
by  the  cathode  plasma.  According  to  the 
Child-Langmuir  current  scaling,  this  gap 
closure  causes  an  increase  in  the  electron 
current  density  by  the  factor  (d  -  ut)“2» 
where  d  is  the  initial  A-K  gap  spacing  and 
u  is  the  closure  velocity.  Another  problem 
intricately  related  to  this  increasing 
current  is  the  maintainance  of  a  flat  output 
voltage  throughout  the  ps  pulse.  Constancy 
of  the  beam  voltage  is  particularly  important 
for  microwave  generators  which  must  satisfy 
resonance  conditions.  The  plasma  closure 
velocity  is  therefore  crucial  in  determining 
diode  impedance  collapse,  voltage  droop,  and 
A-K  gap  shorting. 

MELBA  Design  and  Performance 


The  Michigan  Electron  Long  Beam  Acceler¬ 
ator  (MELBA)  was  designed  (17)  to  generate 
10  kA,  1  MV  electron  beams  into  a  collapsing 
diode  impedance  for  1  ps  pulselengths.  The 
design  utilizes  an  electrical  -voltage  compen¬ 
sation  circuit  first  proposed  by  Abramyan  [9} 
and  later  by  Smith  (10).  This  circuit 


employs  a  Marx  generator  coupled  to  a  ringing 
RLC  circui t . 


Fig.  1.  Simplified  schematic  of  the 

Abramyan  circuit  with  all  spark  gap 
switches  closed  (9,10).  Not  shown 
are  the  Marx/ switch  resistances. 

The  MELBA  design  also  includes  a  small  RC 
filter  circuit  on  the  output.  The  ringing 
atage  (RS)  is  charged  to  initially  deliver 
a  voltage  of  the  opposite  polarity  from  the 
Marx  output.  At  a  later  time,  the  polarity 
of  the  RS  is  in  the  same  sense  as  the  Marx 
output.  Separate  charging  supplies  are 
utilized  for  the  Marx  generator  and  the 
ringing  stage.  This  allows  the  amount  of 
voltage  compensation  to  be  changed  by  choos¬ 
ing  the  relative  charging  voltage  of  the  Marx 
versus  the  ringing  stage.  The  spark  gap 
switches  in  the  Marx  generator  and  the  RS  are 
independently  triggered;  in  these  experiments 
the  two  circuits  were  triggered  simultan¬ 
eously.  MELBA  utilizes  14  capacitors  in  the 
Marx  generator  and  2  capacitors  in  the  RS, 
one  third  the  number  of  capacitors  used  in 
other  microsecond  electron  beam  generators 
12.31 . 

Data  in  Fig.  2  represent  a  -1  MV  out¬ 
put  voltage  from  the  MELBA  machine  into  a  low 
inductance  127  ohm  resistive  load;  also  shown 
is  the  output  of  the  SPICE  transient  circuit 
analysis  code  [IB]  to  illustrate  the  effect 
of  the  ringing,  reverse-charged  stage. 


Fig.  2.  MELBA  output  voltage  into  127  ohm 
resistive  load. 


From  the  above  data  it  can  be  seen  that  the 
experimental  voltage  remaina  flat  to  within 
i  7  %  ( *10  I)  over  pulselengths  of  1  us 
tl.4  /is).  The  tranalant  circuit  code  gives 
excellent  agreement  with  experimental  data. 
It  should  also  be  noted  from  the  computer 
model  that  without  the  ringing  atage,  the 
voltage  would  not  be  flat  but  would  exhibit 
the  RC  decay  characteristic  of  Marx 
rundown  circuits. 

The  generator  was  constructed  with  a 
triggered,  high  pressure  SF.  crowbar  switch 
to  short  circuit  the  generator  output  after 
the  design  pulselength  of  l  at.  The  crowbar 
switch  removed  90*  of  the  generator  voltage 
in  about  40  ns.  Figure  3  gives  the  output 
voltage  and  current  into  a  resistive  load 
without  crowbar  firing. 


load  without  crowbar  firing 
-  (lus/div).  Charging  voltages: 

V  •  ♦  56  kV;  V  -  ♦  54  kV. 

The  data  ox  Fig.  3  demonstrate  the 
effectiveness  of  the  RS  in  generating  flat 
voltage  out  to  1.4  /is.  The  fact  that 
voltage  remains  on  the  load  for  about 
10  /is  into  a  resistive  load  is  important  for 
comparison  with  the  electron  bean  diode 
closure  experiments.  The  resistor  in  series 
with  the  cathode  (Fig.  l)  dissipates  the 
stored  generator  energy  in  the  event  of  a 
crowbar  or  A-K  short  circuit.  This  feature 
permitted  experiments  to  be  performed  without 
firing  the  crowbar  in  order  to  investigate 
plasma  closure  and  shorting  with  multimicro¬ 
second  pulselengths. 

Cathode  -  Anode  Configuration 


in  Fig.  4. 
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Sine*  the  objective  of  these  experiments 
was  to  investigate  high  current  electron 
beams  with  suit  microsecond  pulselengtha,  the 
a^ode-cathode  spacinga  were  relatively  target 
10  cs  for  eost  shots,  although  •  cat  and  6  cat 
were  also  used.  In  order  to  extract  multi- 
kiloampere  currents,  two  types  of  large  area 
field  emission  cathodes  were  eaiployed  in 
these  experiments! 

1)  carbon  brush  cathode,  and 

2)  velvet  cloth  cathode. 

The  brush  cathode  was  constructed 
following  the  design  of  Refs.  4  and  11. 

Carbon  string  was  threaded  through  a  square 
array  of  holes  on  1  cm  centers  in  a  7  cm 
radius  aluminum  plate.  The  outer  radius  of 
the  carbon  brush  array  was  6  cm.  A  smaller 
(2.5  cm  radius)  brush  cathode  was  also 
constructed  with  recessed  CR-39  ion  track 
film  to  detect  ions  crossing  the  A-K  gap. 

The  velvet  cathode  used  in  these 
experiments  was  65  %  cotton  and  35  %  rayon. 

A  6  cm  radius  circle  of  this  cloth  was  epox- 
ied  to  a  7  cm  radius  aluminum  cathode  plate. 

In  both  cases,  the  cathode  plate  was 
threaded  to  the  anodized  cathode  stalk  of 
radius  7  cm.  Anode-cathode  spacinga  were 
adjusted  by  installing  anodised  rings  between 
the  cathode  stalk  and  the  cathode  plate. 

An  extremely  durable  anode  design  was 
required,  since  the  generator  was  typically 
operated  in  a  noncrowbarred  node  in  which 
plasma  shorting  of  the  A-K  gap  occurred.  For 
this  reason  it  was  not  possible  for  foil 
anodes  (or  even  solid  molybdenum  anodes)  to 
survive  the  high  current  discharge  from  the 
capacitor  bank  after  A-K  shorting.  The 
high  current  discharge  was  much  more 
damaging  to  anode  materials  because  high 
current  densities  deposit  energy  at  the 
surface,  removing  material  in  a  vacuum  arc 
mechanism  [19].  It  was  found  that  a  <3.2  mm 
thick)  poco-graphite  anode  produced  the 
least  contamination  of  the  cathode  and 
insulators  for  shorted  A-K  gaps. 

The  anode  plate  was  insulated  from  the 
vacuum  tank  and  up  to  6  current  transf ormers 
were  utilized  to  monitor  total  return  current 
which  passed  through  low  inductance  straps. 
Current  signals  were  summed  in  the  Faraday 
cage  by  ferrite  adders.  Electron  beam 
burrent  was  sampled  by  an  array  of  small 
(1.6  mm)  holes  in  the  graphite  anode  plate. 

A  beam  collector  and  current  transformer 
permitted  measurements  of  apertured  electron 
beam  current.  Energy  discrimination  of 
apertured  electron  beam  current  could  be  per¬ 
formed  by  installing  foil  between  the  carbon 
anode  and  the  beam  collector.  For  studies  of 
electron  beam  uniformity  and  beam  dynamics 
the  electron  beam  collector  was  replaced  by  a 
lucite  plate.  The  inside  of  this  plate 
was  coated  by  an  opaque  layer  of  Aerodag 
graphite  spray >  the  outside  of  this  plate 
was  viewed  by  either  an  open  shutter  or 
streak  camera. 

Electron  Beam  Ceneration  Experiments 


Brush  cathode  turn-on  data  is  presented 
in  the  voltage  and  current  traces  of  Figure  5. 
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Fig.  5.  Cathode  turn-on  data  for  carbon  brush 
with  10  cm  A-K  gap.  (100  ns/div) 
Upper  traces  Voltage  (160  kv/dlv) 
Lower  traces  current  (2.2B  kA/div) 

The  initial  detectable  current  always 
coincided  with  the  knee  (at  about  30  ns) 
during  the  rise  of  generator  voltage.  Rapid 
current  fluctuations  were  measured  for  both 
types  of  cathodes.  The  spiky  nature  of  the 
current  agrees  with  the  NRL  induction  linac 
results  in  which  a  similar  brush  cathode  was 
compared  with  a  thermionic  cathode  [11] . 

The  cause  of  the  current  fluctuations  could 
be  nonsinultaneity  in  the  explosive  emission 
from  the  multiple  tufts  of  carbon  string. 
Voltage  compensation  at  early  times  is 
seen  to  be  excellent,  with  flat  voltage 
(♦  5  t)  over  a  0.7  ^is  pulse. 

Electron  beam  uniformity  was  measured 
by  means  of  the  lucite  plate  diagnostic. 
Figure  6  gives  an  open  shutter  photograph  of 
the  lucite  plate  for  a  MELBA  pulse  which 
crowbarred  at  0.5  ^is. 
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Fig.  6.  Open  shutter  photograph  of  lucite 

plate  for  0.5  its  electron  beam  pulse. 
Brush  cathode  with  10  cm  A-K  gap. 

The  lucite  plate  data  indicates  that  the 
electron  bean  has  good  uniformity  over  the 
full  cathode  diameter.  For  longer  electron 
beam  pulses,  the  large  electrostatic  charge 
deposited  in  the  lucite  caused  arcing 
between  the  apertured  beamlets,  obscuring 
the  data  and  damaging  the  lucite. 

Diagnostic  data  for  two  noncrowbarred 
electron  beam  pulses  from  MELBA  are  presented 
in  Fig.  7. 


Pig.  7  MELBA  electron  beam  pulse  without 

crowbar,  upper  data  represent  Type  A 
shot)  lower  data  set  represents  Type 
B  shot. 

a)  Voltage  (160  kV/divj  0.2/is/div) 

b)  Current 

(5.3  kA/divj  0.2  us/div) 

c)  Voltage  (160  kV/div;0.S  p*/ div) 

d)  Current  (6.2  kA/div:  0.5  ^is/div) 

e)  Hard  X-rays  (0.5  >>s/div) 

f )  Apertured  e-bea*  current 

(20  A/div,0.5  >>s/div) ,  no  foil. 
Parameters*  Brush  cathode  on  7  cm 
radius  plate,  A-K  gap  of  10  cm. 
Charging  voltages: 

V  •  ♦  56  kV  and  V  »  ♦  56  kV 
Two  types  of  current  behavior  were  ob¬ 
served:  Type  A  with  flat  or  decreasing 
current  late  in  the  pulse,  and  Type  B  with 
current  continually  increasing  throughout  the 
pulse.  Por  brush  and  velvet  cathodes,  the 
voltage  flatness  was  limited  by  large  cur¬ 
rent  spikes  and  the  rapid  increase  which 
typically  occurred. 


Large  current  spikes  are  believed  to  be 
cathode  flares  from  individual  tips  and  were 
also  observed  in  lower  current  density  Soviet 
experiments.  The  apertured  e-beam  signal 
exhibited  large  fluctuations  during  enhanced 
current.  These  large  current  fluctuations 
could  be  caused  by  a  streaming  instability. 
(20) 

Diode  Closure  Velocities 


Por  purposes  of  this  research  it  is 
useful  to  define  two  closure  velocities: 

1)  The  effective  closure  velocity  v,  is 
defined  by  the  shorting  time  of  the  A-K  gap. 

2)  Apparent  closure  velocity  u,  is  that  which 

is  defined  in  terms  of  the  Chlld-Langmulr 
current  with  a  closing  A-R  9spi.it  is 
obtained  from  the  slope  of  (P)  1  ,  where 

P  is  the  perveance. 

The  data  of  the  Type  A  shot  in  Pig.  7 
shows  that  the  voltage,  hard  X-rays,  and 
apertured  e-beam  current  persist  for  3.7 
pt  with  a  10  cm  A-K  gapi  this  yields  an 
effective  closure  velocity  of  2.7  em/j>s. 

Figure  •  gives  the  time  history  of  the 
diode  Impedance  and  (P)  ,  for  the  same 

two  shots  given  in  Pig.  7. 
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Figure  I.  Temporal  evolution  of  data  from 
Pig.  7:  a)  Impedance  profile  for 
type  A  shot  (solid  points)  god, type 
B  shots  (open  points),  b)  P 
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i  Not*  that  typ*  A  ahota  gave  a  flat  impedance 

(50  ohms)  after  about  0.1  fit  whereas  type  > 
shots  have  continually  drooping,  impedance. 

In  type  A  shots  the  slop*  of  P  w<  appears  to 
change  at  about  0.1  fit.  This  has  been 
j  interpreted  by  others  as  an  actual  slowing  of 

the  cathode  plasma  velocity  by  *electric 
pressure*  (21]  or  by  the  beat*  Magnetic  field 
(22).  for  a  10  kA  bean  current  the  Magnetic 
pressure  is  Much  less  than  the  kinetic 
pressure  of  typical  cathode  plasma. 

However,  the  relativistic  Child-Langmuir 
analysis  of  the  next  section  shows  that  the 
type  A  current  can  be  Modeled  by  a  constant 
Closure  velocity  at  early  and  late  tines. 

table  1  summarizes  the  apparent  and 
effective  diode  closure  velocities  for  brush 
and  velvet  cathodes  at  different  A-K  gap 
spacing*. 


Table  I 

Summary  of  Pi  ode  Closure  Velocities  (cn/iis) 
BRUSH  VELVET 


A-K  gap (cm)  10 

8 

10 

1 

Apparent 

Type  A  3.3  f.4 

3.3 

Non* 

(3) 

(1) 

Type  B  5.3  s .4 

8 

6.1 

6.2 

(2) 

(1) 

(1) 

(1) 

Effective 

Type  A  3  ♦  .7 

2.1  «.3 

None 

(3) 

(7) 

Type  B  3.5  ♦.« 

2.6 

3  4.6 

4.11.3 

(4) 

(1) 

(4) 

(5) 

(Numbers  in  parentheses  indicate  shots 
averaged) 


It  should  be  noted  that  the  velvet 
cathodes  always  gave  continually  increasing 
current  (type  B)  behavior.  In  spite  of  the 
higher  electric  field  a  number  of  brush 
cathode  shots  with  d  k  ■  S  cm  demonstrated 
lower  closure  velocities  than  the  same 
cathode  with  4|)(  -  10  cm. 

Diode  Current  Data  and  Hodeling 


Our  model  for  current  flow  from  a 
moving  source  plassia  in  the  diode  employs 
a  fully  relativistic  version  of  the  Child- 
Langmuir  current  density  which  can  be 
modified  by  the  inclusion  of  ion  flow  (23] . 
Following  the  technique  of  Poukey,  we  perfora 
a  Romberg  integration  solution  for  the 
current  density  corresponding  to  space  charge 
potentials  in  a  planar  diode  gap.  ror  our 
model  we  have  also  included  closure  of  the  A- 
R  gap  by  the  cathode  plasma  with  velocity  u 
and  ion  flow  of  variable  magnitude  (23-25]. 
Edge  emission  was  Modeled  by  the  technique  of 
Parker  (26). 

Two  types  of  temporal  current  behavior 
were  observed  in  these  experiments.  Figure  9 
compares  type  A  and  type  B  current  traces 
with  theoretical  models. 


Fig. 9  a)  Experimental  current  data  (points) 
from  type  A  shot  compared  with  constant 
closure  velocity  model  (u  -  2.C  cm/iis) 
for  Child-La nqmulr  (CL)  current  and 
bipolar  flow  (BP)  turned  on  linearly 
from  0.45  p*  to  0.7  ^js. 

b)  Experimental  current  data  from  type 
B  shot  conpared  to  Cbild-Langmuir  (CL) 
with  closure  velocity  of  3.1  cm/us  and 
bipolar  (BP)  flow  turned  on  linearly 
from  0.38 /is  to  0.87  fit  with  constant 
closure  velocity. 

c)  Experimental  current  data  from  type 
B  shot  compared  to  Child-Langmuir  (CL) 
with  a  constant  closure  velocity  of  5 
c*/us  and  bipolar  flow  (BP)  turned  on  at 
0.87/is. 


Note  from  rig .  9a  that  the  experimental 
current  aqrees  with  the  Child-Langmuir  model 
from  0  to  0.49  ^aaec  and  from  1.2  >)■  to  1.9 
yt.  However,  between  0.49  fit  and  0.7  fit,  a 
linear  transition  to  bipolar  flow  models  the 
data  very  well.  A  physical  explanation  for 
the  type  A  current  behavior  is  that  at  early 
times  (less  than  0.49  fit)  the  cathode  emits 
the  Child-Langmuir  current.  Between  0.6  yt 
and  0.7  ^is  ions  from  the  anode  or  the 
background  gas  (3  -  8  x  10”3  Torr)  are 
sufficient  to  support  bipolar  flow.  At  later 
times  (greater  than  0.7  jis)  ’ion  starvation* 
causes  the  current  to  return  to  Child- 
Langmuir  flow. 

In  type  B  shots  the  early  time  behavior 
is  siaiilar  to  type  A.  The  type  B  shots, 
however,  exhibit  continually  increasing 
current  throughout  the  pulse.  Figure  9b 
shows  that  even  a  transition  to  bipolar  flow 
(at  constant  closure  velocity)  cannot  model 
this  rapid  current  increase.  An  increased 
closure  velocity  of  9  cm/yt,  modeled  in 
Pig.  9c,  would  provide  a  physical  mechanism 
for  the  rapid  current  increase.  This  agrees 
with  Ref.  27  in  which  the  authors  conclude 
that  for  cases  in  which  a  threshold  energy 
dose  rate  is  exceeded,  the  anode  plastui 
gives  both  a  transition  to  bipolar  flow  and 
an  increased  diode  closure  velocity,  it 
should  be  noted  that  the  energy  density  to 
the  anode  in  our  experiments  is  well  below 
the  required  400-600J/gm  dose  for  graphite 
(27].  The  longer  pulselengths  in  our 
experiments  may  reduce  the  threshold 
energy  for  anode  plasma  production.  For 
microsecond  pulselengths  it  is  also  possible 
that  background  gases  play  a  role  ion 
production. 

The  interpretation  of  enhanced  current 
in  terms  of  anode  plasmas  is  supported  by  the 
observation  that  the  apertured  e-beam  signal 
always  showed  enhanced  fluctuation  levels 
during  the  period  of  increasing  current. 
These  large  fluctuations  could  be  indicative 
of  a  streaming  interaction  between  the  beam 
and  the  anode  plasma  (20].  Etching  of  the 
CR-39  track  film  embedded  in  a  cathode  gave  a 
crater  over  the  entire  exposed  area 
suggesting  a  large  flux  of  low  energy  ions 
backstreamlng  to  the  cathode. 

Some  additional  observations  concern  the 
conditions  under  which  type  A  and  type  B 
currents  were  measured.  Velvet  cathodes 
always  gave  type  B  current.  For  8  cm  A-K 
gaps  the  brush  cathode  was  more  likely  to 
give  type  A  behavior  (7  out  of  6  shots) . 

Brush  cathode  experiments  with  10  cm  A-K 
gaps  showed  that  on  the  first  several  shots 
after  opening  the  diode  to  air,  a  type  A  shot 
was  more  likely >  on  later  shots  a  type  B 
shot  was  more  likely.  Apparently  the  cathode 
current  density,  surface  chemistry  and 
residual  gas  play  a  role  in  microsecond 
plasma  production. 

Transient  Circuit  Codes  and 
Voltage  Compensation 


The  generator  system  performance  was 
calculated  with  the  SPICE  transient  circuit 
analysis  program.  The  computational  model 
Includes  stray  capacitances,  inherent 
inductance  loops,  effective  Marx  shunt 
iiapedance,  and  the  RC  filter  circuit.  The 


values  of  the  various  elements  were  derived 
by  PS1  (17)  from  a  combination  of 
measurements,  experimental  data,  extrapol¬ 
ations,  and  calculations.  The  SPICE  cir¬ 
cuit  code  has  been  run  for  fixed  resistance 
loads  as  well  as  time  varying  impedances. 

For  the  time-varying  diode  impedance,  the 
load  is  modeled  in  a  piecewise  linear  manner 
using  a  nonlinear  circuit  element  (typically 
a  JFET) . 

Current  spikes  and  irreproduclbility  in 
the  diode  current  made  direct  experiMntal 
comparison  of  different  compensating  voltages 
difficult  on  separate  shots. 

Figure  10  compares  the  experimental 
voltage  with  transient  circuit  analysis  code 
results  using  the  experimental  impedance 
profile  for  two  eases: 

1)  Actual  MELBA  circuit  and  charging 
voltages,  and 

2)  MELBA  circuit  but  increased  charging 
voltage  of  ringing  stage. 
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Fig.  10.  MELBA  experimental  voltage  (points) 
for  type  A  shot  compared  to  SPICE  code 
using  experimental  impedance  profile 
with:  Actual  charging  voltages  (dark 
line)  V  •  ♦  56  kV,  V.  .  ♦  56  kV,  and, 
increased  charging  voltage  on  the 
reversed  stage,  V„  »  «56  kV,  V_  »  490  kV 
It  can  be  seen  that  increasing  the  ringing 
reverse  stage  charging  vol'tage  improves  the 
voltage  flatness  between  0.8  us  and  1.4  us. 
However,  the  voltage  deviations  observes  in 
these  experiments  have  been  caused  primarily 
by  current  spikes  and  the  rapid  current 
enhancement  starting  at  0.45  us.  Current 
spikes  may  be  eliminated  by  resistively 
ballasting  the  individual  emitters  on  a 
multi  tip  or  brush  cathode. 

•Supported  by  OMR,  HSF  and  CM  Research  Labs 
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in  Low  Pressure  Media 
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Abstract 

A  Monte  Carlo  Fokker  Planck  code  is  used  to 
investigate  the  scattering  loss  of  particles  from  a 
relativistic  electron  beam  propagating  in  a  low  pressure 
air  channel.  The  fraction  of  particles  transmitted  is 
calculated  and  its  variation  with  relevant  beam  and 
channel  parameters  is  obtained.  The  effect  of  multiple 
scattering  is  examined,  and  a  parameter  reflecting 
scattering  erosion  rate  is  calculated  for  some 
experimentally  observed  parameters  and  found  to  agree 
with  an  analytical  expression  for  force-free  propagation. 


introduction  and  Analysis 

Some  of  the  most  comprehensive  results  on  beam 
erosion  in  the  Ion-Focused  Regime  (IFR)  have  recently 
emerged  from  the  ETA  experiment  at  Livermore.  '  In  air 
this  regime  for  stable  propagation  corresponds  to  a 
pressure  window  of  about  0.1-0. 2  torr.  At  0.08  torr  beam 
erosion  primarily  at  the  head  of  the  beam  was  observed  and 
found  to  correspond  to  0.16  m/m.  At  these  low  pressures 
it  is  expected  that  most  of  the  erosion  will  be  due  to 
ohmic  heating;  however  it  is  also  important  that  some 
assessment  of  the  role  of  scattering  in  erosion  be  made 
particularly  as  the  parameters  of  the  beam  or  the  channel 
vary.  In  this  paper  we  examine  the  loss  of  particles  from 
the  beam  as  it  transverses  a  preformed  channel  and  attempt 
to  deduce  quantities  which  may  have  a  direct  bearing  on 
erosion.  The  investigation,  is  carried  out  with  the  aid  of 
a  Monte  Carlo  Fokker  Planck  code  which  allows  us  to  follow 
a  given  sample  of  non-interacting  beam  particles  in  full 
configuration  and  velocity  space  as  the  beam  propagates  in 
a  cylindr ically  symmetric  gaseous  channel.  It  is  assumed 
that  the  channel  medium  is  partially  ionized  so  that 
interaction  of  the  beam  particles  with  charged  and  neutral 
target  particles  are  taken  into  account.  The  probability 
of  neutral  particle  interaction  is  selected  by  defining  an 
appropriate  mean  free  time 


where  17  is  the  beam  electron  velocity,  n  is  the  neutral 
density  and  0“  is  an  appropriate  interaction  cross 
section.  Although  other  inelastic  events  should  generally 
be  considered  we  limit  our  study  to  ionization  and  the 
interaction  probability  is  given  by  dt/tf  when  dt  is  a 
very  small  code  time  step.  Given  an  appropriate 
differential  scattering  cross  section  a  random  scattering 
angle  £?  can  be  selected  from  the  distribution 
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where  X  is  -a  random  number.  For  relativistic  electrons 
propagating  in  a  low  temperature  medium  it  is  expected 
that  the  distribution  be  highly  peaked  in  the  forward 


direction  and  using  the  results  from  hydrogen 
appropriately  scaled  for  higher  z  elements  we  can  write 
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where  in  CGS  units 


[f  V  '  1 

In  the  above  expression  m  is  the  electron  mass,  k  is  the 
Debroghi  wave  number,  h  is  Plancks  constant,  all. 4  a  z~ 

7  ,  a  is  the  Bohr  radius  and  &  •  is  the  minimum 
scattering  angle  given  by 


mmnmiiii 


(3) 


where  Y  is  the  relativistic  parameter  and  c  is  the  speed 
of  light.  Although  small  angle  Coulomb  scattering  is 
included  in  the  code  by  an  implicit  Fokker  Planck  method 
the  probability  for  such  scattering  event  in  the  code  time 
step  dt  is  simply  dt/’C-  where  7,  i6  the  well  known 
deflection  time.  For  plasma  densities  on  the  order  of 
1(T  -lO^cm”'3  which  are  appropriate  for  the  case.at  hand, 
and  for  tT  *  c  it  is  easily  seen  that  Td  s  10*"A  seconds 
so  that  at  pulse  lengths  on  the  order  of  nanoseconds  it  is 
clear  that  scattering  from  charged  particles  can  be 
neglected.  These  considerations  are  consistent  with  some 
experimental  observations. 1,91 

In  order  to  crudely  model  the  propagation  of  a 
relativistic  electron  beam  in  either  an  IFR  or  a  force- 
free  propagation  mode  we  have  incorporated  the  self 
magnetic  field  in  a  phenomenological  way  using  Ampere's 
law  for  an  infinite,  uniform  beam.  We  have  neglected  the 
radial  electric  field  by  taking  a  charge  neutralizing  ion 
density  (n-)  to  be  equal  to  the  beam  density  (n.)  so  that 
the  chargeAneutralization  parameter  f  -  1.  Since  the  code 
calculates  the  number  of  particles  that  traverse  the 
channel  the  fraction  of  the  beam  that  is  transmitted  would 
then  represent  an  upper  limit  for  pulse  lengths  7"  that 
are  shorter  than  XCN#  the  charge  neutralization  time, 
and  a  good  estimate  when  *ZL  -  TCN  since  the  physics  of 
instabilities  that  might  occur  are^not  included.  Under 
these  conditions  the  description  of  the  scattering  process 
might  be  deemed  appropriate 2to  that  taking  place  in  a 
force  neutralized  f  >  1/ y  ,  charge  neutralized,  f  “  1 
or  force  free  f  ~  I,  f  “  1  regimes.  With  propagation  _ 
having  been  obslrved  in  pressure  ranges  of  0.1-5  torr14'1” 
it  should  be  noted  that  at  lower  pressures  there  are 
usually  insufficient  neutrals  to  provide  an  adequate  ion 
density  for  force  neutralization  and  at  higher  pressures 
multiple  scattering  could  cause  rapid  dispersal  of  the 
beam.  The  latter  effect  is,  however,  incorporated  in  the 
code  and  its  influence  will  be  discussed  later. 

In  assessing  the  utility  of  the  Monte  Carlo-Fokker 
Planck  code  in  the  beam  propagation  study  it  should  be 
kept  in  mind  that  particle  scattering  can  be  accurately 
described  and  that  no  assumptions  need  be  made  concerning 
the  plasma  conductivity  since  ionization  can  be  calculated 
self-consistently.  Particle  and  energy  loss  mechanisms 
can  also  be  readily  included  so  long  as  the  appropriate 
mean  free  paths  for  these  processes  can  be  selected. 
Because  of  its  three  dimensional  simulation  capability  the 
often-used  "paraxial"  approximations  and  others  need  not 
be  made.  Because  of  its  limitation  in  the  number  of 
particles  that  can  be  simulated  it  is  clear  that  the 


treatment  of  self  fields  and  the  modelling  of 
instabilities ,  though  possible,  might  be  quite  difficult 


Boms.  Preliminary  Results 


A  relevant  factor  to  beam  erosion  by  scattering  is 
whether  the  beam  particle  undergoes  single  or  multiple 
scattering  events  as  it  transverses  the  channel.  If  we 
denote  by  aTn  the  transmitted  fraction  of  the  beam,  n  the 
number  of  times  scattering  off  of  a  neutral  particle  m 
occurs,  n  the  number  of  electrons  used  to  simulate  the 
beam,  r  the  radius  of  the  channel  and  s  =  n  /n  then  the 
variation  of  these  parameters  with  r  for  a  DeafiLradius 
r.  =  1.5  cm,  total  cross  section  cr  £  *  =  3xlO-ABcm  » 
neutral  density  n  *  2xlO,i:>cm-J  and  n°  =  5000  is  shown  in 
Table  1.  n 
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0.007 

7338 

1.47 
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'•EH29H 

9258 

1.85 

0.03 

12183 

2.4 

15 

0.079 

15003 

3.0 

21.2 

0.14 

17800 

3.56 

We  observe  that  the  transmitted  fraction  is  linearly 
proportional  to  r*  which  is  a  direct  consequence  of  the 
geometric  effect.  cAnother  geometric  consequence  is  the 
linear  scaling  of  n  and  s  with  r  .  Since  the  probability 
of  neutral-interaction  is  proportional  to  the  path  length 
traversed,  and  the  path  length  increases  with  r  then  the 
probability  of  multiple  scattering  should  increase  with 
r  .  The  variation  of  these  parameters  with  the  neutral 
density  is  shown  in  Table  2  for  a  channel  radius  r  =  21.2 
cm  and  beam  radius  and  total  scattering  cross  section  as 
before. 
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293 

0.06 

0.95 

2xl014 

2789 

0.66 

0.79 

2xl015 

18000 

3.6 

0.14 

2xl016 

68263 

13.65 

0.002 

The  transition  from  single  to  multiple  scattering#  s  s  1# 
is  clearly  noted  in  the  table.  Although  the  distribution 
for  single  scattering  is  highly  peaked  in  the  forward 
direction#  beam  particles  may  still  gain  sufficiently 
large  transverse  directed  energy  to  cause  substantial  beam 
breakup. 

The  code  was  also  utilized  to  estimate  an  "erosion 
rate*  due  to  scattering.  It  yields  the  fraction  of  the 
beam  lost  per  unit  length  travelled.  When  this  fraction 
is  multiplied  by  the  length  of  the  beam  a  loss  rate  due  to 
scattering  is  obtained.  This  is  given  by 
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where  m/*m  is  the  fraction  of  particles  lost  per  meter# 
and  LTTd  *  I*  is  the  length  of  the  pulse  in  meters. 
Equation  (4)pis  based  on  the  assumption  that  for  constant 
beam  density  and  constant  loss  rate  of  particles  the  loss 
rate  of  particles  in  a  rigid  beam  is  equivalent  to  the 
loss  rate  of  beam  length.  When  the  parameters  of  the 
force  free  propagation  experiment  of  Ref.  4  are  used#  Eq. 
(4)  yields  results  comparable  to  those  predicted  by  the 
analytical  expression  given  by  Leelo;  that  is  based  on  the 
SNEEZE  code. 
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Gas  Hydrodynamics  of  Channel*  formed  bjr 
Rapid  Heatlnq  by  Lasers  and  Later  Guided 
PlscKarges, *  L.D.  HORTbN  and  R.M. 
GILGENBACH,  Department  of  Nuclear 
Engineering,  The  University  of  Michigan, 
Ann  Arbor,  Ml,  48109. 

The  hydrodynamic  evolution  of 
reduced  density  channels  produced  by  CO, 
laser  induced  gas  breakdown  and  by  laser 
guided  discharges  have  been  studied.  In 
the  laser  induced  breakdown  experiments, 
a  high  power  (0.30  GW)  TEA  CO,  laser  is 
brought  to  a  quasi-line  focus  over  14  cm 
in  argon  by  a  compound  reflecting  axicon 
with  a  w-shaped  cross  section  (waxicon) . 
Electrical  discharges  of  30  kV  are 
initiated  and  guided  by  CO,  laser 
induced  air  breakdown  focused  by  a 
conventional  mirror  with  2.5  m  focal 
length. 

Experimental  data  are  obtained 
using  schlieren  photography  and 
col  linear  holographic  interferometry. 
Shock  positions  and  directly  measured 
(non-Abel  inverted)  channel  density 
profile  data  are  compared  to  a  one 
dimensional  hydrodynamics  code  in  order 
to  estimate  fractions  of  available 
energy  absorbed  into  channel  heating. 
This  one  dimensional  code  is  valid  for 
times  short  compared  to  the  onset  time 
for  convective  mixing  which  eventually 
cools  and  breaks  up  the  reduced  density 
channel. 

For  waxicon  focused  laser  induced 
gas  breakdown  channels  the  measured 
shock  wave  velocity  ranges  from  0.4  to 
0.5  mmlfUS.  Comparison  to  the 
computational  model  indicates  that  25- 
301  of  the  laser  energy  is  absorbed  into 
heating  the  channel. 

For  laser  guided  discharges  the 
measured  density  at  the  center  of  the 
channel  is  0.025  kg/m  .  This  is  in 
excellent  agreement  with  the  one 
dimensional  code  for  times  small 
compared  to  the  time  scale  for 
convective  mixing.  However,  at  later 
times  the  code  predicts  narrower 
channels  with  lower  on-axis  densities 
than  are  measured. 

•This -research  is  supported  by  the 
Office  of  Naval  Research,  National 
Science  Foundation  Grant  ECS-8309682, 
and  a  Presidential  Young  Investigator 
Award  funded  by  General  Motors  Research 
Labs  and  the  National  Science 
Foundation. 
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Spatially  and  Temporally  Resolved  Diagnostic* 
for  Microsecond,  Intense  Electron  Beama  * 

R.  M.  GILGENBACH,  M.  BRAKE,  L.  D.  HORTON, 

S.  B I  DWELL,  R.  F.  LUCEY,  L.  SMUTEK,  and 

J.  E.  TUCKER,  Nuclear  Engineering  Department, 
University  of  Michigan,  Ann  Arbor,  Ml  48109 
Experiments  are  underway  to  investigate  new 
diagnostics  for  electron  beams  in  vacuum  and 
in  a  plasma  background.  Measured  parameters 
include  temporally  resolved  beam  current 
profile  and  beam  emittance.  These  character¬ 
izations  are  being  performed  during  electron 
beam  diode  closure  experiments  (1)  and  beam- 
plasma  interaction  experiments  with  either  of 
two  long-pulse  accelerators:  MELBA  (Michigan 
Electron  Long  Beam  Accelerator) : 

Voltage  •  -  1  MV,  Current  *  10  kA,  at 
Pulselength  ■  0.1  to  lxs  (1.4  MS  )  for 
voltage  flat  to  within  t  7%  ( ±  10%)  . 

The  second  accelerator  is  a  long-pulse 
Febetron  with  parameters:  Voltage  ■  -  0.5  MV, 
Current*  1  kA,  and  Pulselength  *  0.3  s. 

Two  different  configurations  have  been 
developed  which  use  Cerenkov  radiation  to 
detect  electron  beam  current  profiles  as  a 
function  of  time.  The  first  uses  Cerenkov 
emission  by  electrons  which  impinge  axially 
on  a  single  fiberoptic  lightguide  enclosed 
in  a  lucite  tube.  Plasma  light  is  blocked 
by  graphite  spray  or  thin  foil  covering  the 
end  of  the  optical  fiber.  This  diagnostic 
has  the  the  following  advantages: 

1)  the  threshold  energy  for  Cerenkov  emission 
effectively  discriminates  between  high  energy 
beam  electrons  and  low  energy  (3-5  eV) 
plasma  electrons.  2)  The  small,  nonconducting 
probe  introduces  a  minimal  perturbation  into 
the  beam-plasma  system,  3)  Excellent  signal 
to  noise  ratio  is  obtained  because  the 
fiberoptic  signal  is  directly  transmitted 
to  a  photomultiplier  tube  in  the  Faraday  cage. 
4)  Quantitative  data  is  obtained  directly. 

The  second  Cerenkov  diagnostic  employs 
an  array  of  small  apertures  in  the  anode 
plate.  The  apertured  beamlet  array  strikes  a 
clear  Cerenkov  plate  located  several  cm  behind 
the  anode.  This  plate  is  covered  by  a  thin 
layer  of  graphite  to  pass  electrons  but  not 
diode  light.  The  back  of  the  Cerenkov  plate 
can  be  viewed  by  open  shutter  or  streak 
photography.  This  diagnostic  thus  has  the 
potential  for  temporally  resolved  measurements 
of  both  beam  current  density  and  emittance. 

Time  resolved  optical  emission  spectro- 
copy  is  used  to  determine  plasma  composition, 
degree  of  ionization,  and  equilibrium  con¬ 
ditions  of  long  pulse  REB  produced  plasmas 
at  pressures  above  0.1  Torr.  A  gated  inten¬ 
sified  diode  array  detector  is  employed  with 
a  275mm,  f/3.8  spectrograph. 

*  This  research  is  supported  by  the  Office 
of  Naval  Research  and  a  Presidential  Young 
Investigator  Award  Funded  by  the  National 
Science  Foundation  and  General  Motors 
Research  Laboratories. 
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Microsecond  electron  beam  experiments  are 
being  conducted  on  MELBA  (Michigan  Electron 
Long  Beam  Accelerator).  This  new  Long  Pulse 
accelerator  operates  with  parameters: 
Pulselength  =  0.1  to  1.4  microseconds 
Voltage  =  -1  MV  flat  to  <  +  10  %  at 
Current  =  10  kA. 

Unique  circuitry  permits  electrical 
compensation  for  diode  impedance  droop. 

Initial  interaction  experiments  used  Febetron 
beams  with  0.3  microsecond  pulsewidths  at 
voltages  <  500  kV  and  currents  <  1.4  kA. 
Interaction  data  with  argon  and  nitrogen  at 
0.5  Torr  to  40  Torr  show  induced  plasma 
currents  which  persist  for  >  1  microsecond 
after  the  beam  pulse  ends.  Argon  gives  longer 
decay  times  than  nitrogen.  New  Industrial 
processing  applications  of  long-pulse  beams 
are  being  explored. 
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ABSTRACT 

Reduced  density  channels  produced  by  laser  induced  gas 
breakdown  and  by  laser  guided  discharges  are  modelled  by 
means  of  a  one  dimensional  gas  hydrodynamics  code  for  times 
earlier  than  the  onset  of  convective  mixing.  Comparison  of  laser 
induced  gas  breakdown  experiments  with  the  code  indicate  that 
the  conversion  of  laser  energy  to  gas  heating  is  limited  to  25- 
30%.  Col linear  holographic  measurements  of  channel  density 
profiles  produced  by  laser  guided  discharges  are  found  to  be  in 
good  agreement  with  code  predictions  for  times  less  than  50  /ns. 
This  agrees  with  other  experiments  which  have  indicated  onset 
times  for  convective  mixing  in  the  range  50-70 after  the  start 
of  the  discharge. 


Introduction 


Extensive  experiments  at  the  Naval  Research  Labs  (NRL)  have 

studied  channels  formed  by  the  rapid  deposition  of  energy  in 

atmospheric  pressure  gas.1-6  Laser  guided  discharges  and  the 

resulting  reduced  density  channels  have  applications  in  electron 

and  light  ion  beam  transport3' 7-9 ,  in  lightning  modelling10,  and 
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in  materials  processing  '  .  Here  we  report  results  for  the 

nonturbulent  phase  of  reduced  density  channels  formed  by  laser 

induced  gas  breakdown13  and  by  laser  guided  electrical 
14 

discharges  .  Data  is  analysed  and  interpreted  using  a  one 
dimensional,  cylindrically  symmetric  gas  hydrodynamics  code. 

This  code  uses  the  standard  conservation  equations  in  Lagrangian 
coordinates  with  the  addition  of  an  artificial  viscosity15  to 
treat  shock  waves. 

Computations  by  Picone  et  al4'16^17  at  NRL  have  shown  that 
vorticity  leads  to  turbulent  convective  mixing10  which  is  the 
dominant  decay  mechanism  for  reduced  density  channels.  Since  our 
one  dimensional  hydrodynamic  model  cannot  include  vorticity,  we 
model  only  the  pre-mixing  phase  as  determined  from  schlieren  and 
holographic  data.  Col linear  holographic  interferometry  is  used 
in  our  experiments  to  directly  measure  the  density  profiles 
existing  within  the  channel  for  comparison  with  the  code. 

I  Experimental  Configurations 

Two  different  techniques  are  employed  for  rapidly  depositing 
energy  in  a  gas  channel: 

1)  C02  laser  induced  gas  breakdown  (Fig.  la) 

2)  laser  guided  discharges  (Fig.  lb) 


The  COj  laser  gas  breakdown  experiment  uses  a  waxicon  to 
produce  a  quasi-line  focus13  as  depicted  in  Fig.  la.  A  waxicon 
is  a  focusing  mirror  with  a  W  shaped  cross  section.  Argon  is 
used  as  the  gas  for  the  waxicon  experiments  reported  here.  The 
schlieren  diagnostic  system  employs  a  pinhole  aperture  to  obtain 
two  dimensional  spatial  resolution. 

The  laser  guided  discharge  experiments  employ  a 
configuration  (Fig.  lb)  which  permits  the  interferometric  data  to 
be  taken  in  an  orientation  col linear  to  the  discharge.  A  high 
power  (0.30  GW)  TEA  C02  laser  is  focused  by  a  2.5  m  focal  length 
gold  mirror.  High  voltage  electrodes  are  located  across  the 
length  of  the  resulting  breakdown  beads.  A  0.04  /i?  capacitor 
charged  to  30  kv  is  discharged  across  the  electrode  gap  5 
after  the  C02  laser  pulse.  Discharge  current  is  a  damped 
sinusoid  with  peak  current  of  about  10  kA,  period  of  1.4 /as,  and 
1/e  decay  time  of  3.5  The  uniformity  and  decay  time  of 

channels  generated  by  this  experiment  have  been  demonstrated  to 
closely  resemble  the  longer  channels  at  NRL3.  Two  pulse 
holographic  interferometry  is  used  to  record  the  data.  Between 
pulses  the  precision  mirror  is  rotated  slightly  to  produce 
secondary  interference  fringes  on  the  film.  Deflection  of  these 
fringes  by  the  object  gives  a  measure  of  the  gas  density.  Due  to 
the  length  of  the  object  when  the  channel  is  viewed  in  this 
orientation  and  the  resulting  large  fringe  deflections,  helium, 
which  has  a  low  index  of  refraction,  is  used  as  the  ambient  gas. 
Note  that  Abel  inversion  of  holographic  data  is  not  required  in 
the  col linear  arrangement  since  the  data  directly  yields  the 


density  profile  in  the  interior  of  the  channel 
II  Numerical  Model 


The  gas  hydrodynamics  code  uses  the  standard  conservation  of 
mass,  momentum,  and  energy  equations  in  Lagrangian  coordinates 


and  cylindrical  symmetry. 


dp  JjL 

dn  _  _  JL 
dt  P  dr 


d_2_  =  - ( ^ p  +  ‘r 

d-e  v  ° 

<dd 

Here  p  is  the  gas  pressure,  ^  is  the  gas  density,  u  is  the  local 
velocity,  and  ^ is  the  ratio  of  specific  heats.  Note  that  q  is 
an  artificial  viscosity,  originally  proposed  by  VonNeumann  and 
Richtmyer15,  which  is  taken  to  be  quadratic  in  the  time 
derivative  of  the  density  so  that  discontinuites  at  shock  fronts 
are  smoothed  over  a  few  mesh  spacings.  Using  the  continuity 
equation  the  artificial  viscosity  takes  the  form 


Here  c  is  a  constant  varied  as  a  parameter  to  ensure  smooth 
behavior  around  the  shock.  For  the  results  given  in  this  paper 
c  varies  in  the  range  2  to  5.  An  explanation  of  the  logic 
leading  to  this  form  of  q  is  given  in  Ref.  19.  However,  q  is, 
in  fact,  artificial  and  other  variations  have  been  employed. 

The  difficulty  in  modelling  these  experiments  lies  in  the 
driving  term.  Standard  blast  wave  theory  assumes  a  source  which 
is  a  delta  function  in  time  and  simply  sets  initial  conditions 
corresponding  to  an  appropriate  overpressure.  This  approach  is 
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acceptable  for  our  laser  induced  gas  breakdown  experiments  but 
the  laser  guided  discharge  experiment  has  an  energy  deposition 
time  which  is  comparable  to  hydrodynamic  time  scales.  The  time 
dependence  of  the  energy  deposition  for  the  laser  guided 
discharge  is  found  by  using  the  current  histories  and  by  assuming 

that  the  channel  resistance  is  constant  so  that  E=Rll  dt.  The 

■'o 

actual  current  traces  closely  follow  the  form  of  a  damped  sine 
wave  which  agrees  with  the  assumption  of  a  constant  resistance. 
Since  the  total  absorbed  energy  is  the  parameter  being  varied  in 
the  code,  the  absolute  value  of  the  channel  resistance  is  not 
required. 

The  spatial  dependence  of  the  energy  deposition  is  not 
straightforward.  For  the  waxicon  experiment  the  minimum 
breakdown  spot  size  is  limited  by  machining  irregularities  in  the 
mirror  and  is  estimated  by  schlieren  photographs  taken  just  at 
the  end  of  the  CC>2  laser  pulse  to  be  0.2  mm  in  diameter. 

Doubling  this  spot  size  has  only  a  small  effect  on  the  code 
output  and  does  not  change  the  estimates  of  absorbed  energy 
fraction  reported  below. 

The  laser  guided  discharge  energy  deposition  radius  is  a 
function  of  time  due  to  channel  expansion  during  the  discharge. 
The  channel  radius  is  estimated  using  the  data  in  Fig.  2  for  an 
air  channel.  Detailed  transverse  schlieren  photographs  for 
atmospheric  pressure  helium  are  difficult  to  obtain  at  the  ruby 
laser  probe  wavelength  (694.3  nm)  but  schlieren  data  does 
indicate  channels  which  are  about  25%  wider  in  helium  than  in 
argon,  N2,  or  air.  This  factor  is  incorporated  into  the  code's 
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estimate  of  radius  of  energy  deposition.  While  it  is 

conceptually  preferable  to  let  the  code  calculate  the  channel 

radius  from  the  density  data  of  the  previous  time  step.  Fig.  2 

shows  that  the  channel  expands  quickly  to  the  laser  breakdown 

bead  radius.  The  initial  breakdown  phase  is  not  modelled  in  the 

code  for  laser  guided  discharges  and  is  clearly  not  one 

dimensional.  Therefore,  an  external  input  of  early  time  channel 

radii  is  necessary.  Note  that  channel  data  from  over  a  thousand 

shots  show  reproducible  straight,  approximately  one  dimensional 

channnels  at  times  greater  than  about  10  s  after  the  beginning 

of  the  discharge  (see  the  photographs  in  Fig.  2  for  31  and 

51  yws) .  This  indicates  initial  laser  asymmetries  are  not 

important  in  the  nonconvective  region  of  laser  guided  discharge 

3_6  1  ft 

channel  evolution  ' 

III  Results  and  Discussion 

A  comparison  of  experimental  and  code  shock  positions  for 

the  waxicon  breakdown  experiment  is  shown  in  Fig.  3.  The 

available  energy  for  the  experimental  points  is  calculated  by 

using  laser  calorimetry  to  normalize  the  power  density  on  axis  as 

13 

given  by  a  ray  tracing  code  .  As  described  in  reference  13  the 
waxicon  mirror  produces  one  region  on  axis  where  singly  and 
doubly  reflected  rays  cross  and  another  region  with  only  doubly 
reflected  rays  crossing  the  axis.  Both  of  these  regions 
experience  approximately  constant  energy  depositions  and 
correspond  to  approximately  straight  shock  fronts.  Thus,  two 
data  points  are  available  for  each  incident  laser  energy.  Shock 
positions  are  very  reproducible  for  a  given  laser  energy  setting. 
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The  estimate  of  available  energy  is  undoubtedly  better  for  the 

longer  region  of  only  doubly  reflected  rays  but  both  energy 

estimates  have  errors  of  less  than  10%.  Code  results  for  various 

absorbed  energies  are  superimposed  on  the  experimental  results  in 

Fig.  3.  Comparison  between  the  data  and  the  code  indicates  that 

only  25-30%  of  the  available  laser  energy  is  going  into  heating 

20 

the  gas.  Several  potential  energy  loss  mechanisms  exist 
Ionisation  of  the  gas  followed  by  recombination  radiation  which 
escapes  the  channel  is  one  possibility  but,  for  our  short  laser 
pulse  length  and  small  focal  spot  size  this  is  expected  to  be 
negligible.  Inaccuracies  in  the  focusing  mirror,  producing  off 
axis  rays  is  another  potential  mechanism.  The  focal  radius  as 
predicted  by  early  time  schlieren  photographs  suggests  that  this 
also  is  a  small  contribution  to  the  large  energy  fraction  not 
heating  the  gas.  The  most  likely  mechanism  for  major  energy 
loss  is  supercritical  plasma  reflection.  Once  the  breakdown 
channel  is  completely  singly  ionized  by  the  laser  the  resulting 
plasma  is  above  critical  density  for  10.6 y<m  propagation  and  the 
remaining  light  is  reflected  rather  than  being  absorbed.  Only 
nonlinear  processes  can  heat  the  channel  once  it  has  reached 
critical  density  which  explains  the  poor  coupling  efficiency  seen 
in  Fig.  3. 

The  experimental  data  for  the  laser  guided  discharges  are 
obtained  using  collinear  holographic  interferograms  such  as  the 
one  in  Fig.  4.  Relative  densities  are  calculated  by  measuring 
fringe  deflections  from  the  reference  at  the  edge  of  the  channel. 
This  reference  density  was  originally  estimated  to  be  ambient  but 
code  results  indicate  a  slight  density  reduction  behind  the  shock 
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especially  for  times  less  than  50 /is.  Thus  the  experimental  data 

for  Fig.  5  is  normalized  to  the  code  prediction  of  the  density 

just  outside  the  channel.  The  code  best  reproduces  the 

experimental  results  for  an  absorbed  energy  of  100  J/m.  The 

available  energy  can  be  roughly  estimated  to  be  110  J/m  by  noting 

that  the  final  capacitor  voltage  after  the  discharge  is  16  kV  and 

by  assuming  about  half  of  the  discharged  capacitor  energy  is 

dissipated  in  the  spark  gap  switch.  Agreement  for  channel 

density  profiles  is  very  good  at  30  /.is  but  for  times  later  than 

50 yMS  convective  mixing  has  become  important  and  the  code 

predicts  deeper,  narrower  reduced  density  channels  than  are 

actually  measured.  This  onset  time  for  convective  mixing  agrees 

quite  well  with  data  from  transverse  schlieren  photographs  which 

begin  to  show  enhanced  line  integrated  turbulence  at  about  70yMS 

18 

after  the  start  of  the  discharge  . 

Fig.  6  shows  the  pressure  profile  corresponding  to  the 
reduced  density  channel  of  Fig.  5(a).  Note  that  by  30,ms  the 
channel  has  already  expanded  to  pressure  equilibrium  and  has 
essentially  finished  its  evolution  in  our  convectionless, 
conductionless  equations. 
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Figure  Captions 

Fig.  1.  (a)  Experimental  configuration  for  the  waxicon  focused 

C02  laser  induced  gas  breakdown  experiment.  (b)  Experimental 
configuration  for  col  linear  holographic  interferometry  of  laser 
guided  discharges. 

Fig.  2.  Schlieren  photographs  at  different  times  during  the 
evolution  of  a  laser  guided  discharge.  Upper  times  are  relative 
to  the  guiding  laser  pulse,  lower  times  to  the  start  of  the 
electrical  discharge. 

Fig.  3.  Comparison  of  theory  and  experiment  for  shock  positions 
in  waxicon  focused  C02  laser  induced  gas  breakdowns.  Points  are 
experimental  data  and  curves  are  code  output.  =  53  J/m 

available  energy,  (Q)  =  91  J/m,  (O) =  105  J/m,  (0)*  180  J/m;  dot- 
dashed  curve®  10  J/m  absorbed  energy,  long  dashed  curve®  30  J/m, 
short  dashed  curve®  50  J/m,  solid  curve®  70  J/m. 

Fig.  4.  Collinear  holographic  interferogram  for  a  6  cm  helium 
laser  guided  discharge  channel  31  y+s  after  discharge  initiation. 

Fig.  5.  Density  profiles  for  6  cm  helium  laser  guided 
discharges.  Points  are  experimental  data  and  curves  are  one 
dimensional  hydrodynamics  code  results.  (a)  31yns  after 
discharge,  (b)  51  (c)  71  y*s. 

Fig.  6.  Pressure  profile  for  a  6  cm  laser  guided  discharge 
channel  in  helium  31 after  discharge  initiation. 
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COMPARISON  OF  THE  L/R  DECAY  TIMES  FOR 
THE  TWO  6ASES  INDICATE  THAT  THE 
POST-BEAM  RESIDUAL  RESISTIVITY  OF 
NITROGEN  PLASMA  CHANNELS  IS  AT  LEAST  A 


